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Abstract 

Recently, detailed calculations of the excitonic insulator phase model adapted to the case of 1 T-TiSe2 have been presented. 
Through the spectral function theoretical photoemission intensity maps can be generated which are in very good agreement with 
experiment [Phys. Rev. Lett. 99, (2007) 146403]. In this model, excitons condensate in a BCS-like manner and give rise to a charge 
density wave, characterized by an order parameter. Here, we assume an analytical form of the order parameter, allowing to perform 
temperature dependent calculations. The influence of this order parameter on the electronic spectral function, to be observed in 
photoemission spectra, is discussed. The resulting chemical potential shift and an estimation of the resistivity are also shown. 

Key words: photoemission, strongly correlated electrons, calculated intensity maps, exciton condensate 
PACS: 71.45.Lr, 71.27. +a, 79.60.Bm, 71.35.Lk 
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1. Introduction 

The transition metal dichalcogenide 1 T-TiSe2 is a lay- 
ered compound exhibiting a commensurate (2x2x2) charge 
density wave (CDW) at low temperature [Tj. In transport 
measurements, its signature is a strong anomaly in the 
resistivity, peaking slightly below the transition temper- 
ature T c = 200K. The CDW is accompanied by a pe- 
riodic lattice distortion involving small atomic displace- 
ments. At the same time, a zone boundary phonon softens, 
being a central ingredient to conventional theories of struc- 
tural transitions [2J. High resolution angle-resolved photoe- 
mission spectroscopy (ARPES) measurements brought a 
deeper insight into the electronic properties of 1 T-TiSe2 by 
clearly revealing the appearance of new backfolded bands, 
characterisitc of a new periodicity in the system j3|4|5| . 

The origin of the CDW can hardly be explained by the 
usual nesting mechanism [6j, because there are notably no 
large parallel Fermi surface portions [7J. Currently, the best 
candidates are a band Jahn- Teller effect [8J and the exci- 
tonic insulator phase. Recently, the latter scenario has been 
strongly supported by comparison of ARPES data with 
theoretical photoemission intensity maps [9]. These calcu- 
lations are based on the excitonic insulator phase model 
|10)llj . which has been adapted to the particular case of 
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1 r-TiSe2p2]. The excitonic insulator phase may occur in a 
semimetallic or semiconducting system exhibiting a small 
(negative respectively positive) gap. Indeed, for a low car- 
rier density, the Coulomb interaction is weakly screened, al- 
lowing therefore bound states of holes and electrons, called 
excitons, to build up in the system. If the binding energy Eb 
of such pairs is larger than the gap Eq , the energy to create 
an exciton becomes negative, so that the ground state of the 
normal phase becomes unstable with respect to the sponta- 
neous formation of excitons. At low temperature, these ex- 
citons may condense into a macroscopic coherent state in a 
manner similar to Cooper pairs in conventional BCS super- 
conductors. Exciton condensation may lead to the forma- 
tion of CDW of purely electronic origin (not initiated by a 
lattice distortion), characterized by an order parameter A. 
To our knowledge, lT-TiSe2 is the only presently known 
candidate for a low temperature phase transition to the ex- 
citonic insulator state without the influence of any external 
parameters other than temperature. Indeed, as pressure is 
increased above 6 kbar on TmSeo.45Teo.55 ( controlling the 
gap size and thus the energy necessary to create excitons), 
a transition to an insulating phase happens, whose origin 
can also be explained with exciton condensation [13] . In this 
context Bronold and Fehske proposed an effective model 
for calculating the phase boundary of a pressure-induced 
excitonic insulator, in the spirit of a crossover from a Bose- 
Einstein to a BCS condensate 14 1. 



Here, we study the temperature dependence of the exci- 
tonic insulator phase. This is achieved by choosing a simple 
analytical form for its order parameter and inserting it into 
the results of our previous calculations [12] . 

2. Results and Discussion 

The electronic structure of 1 T-TiSe2 near the Fermi en- 
ergy Ep is composed of three Se4p-derived valence bands 
at the T point (center of Brillouin zone) and three Ti3<i- 
derived conduction bands distributed among the three 
symmetry equivalent L points (zone boundary). There is 
a slight overlap of ~ 70 meV, such that lT-TiSe2 has 
a semimetallic character. This issue is still controversial 
among the ARPES community, but a recent infrared study 
confirmed the semimetallicity [15] . In an ionic picture, 
Ti [Ar]3<i 2 4s 2 gives all its valence electrons to the two 
neighbouring Se [Ar]3d 10 4«s 2 4j9 4 , leaving a system with an 
empty d° shell. In our model, we consider only the topmost 
valence band (the other two do not cross the Fermi energy 
and play a minor role) and the three conduction bands. 
The valence band and the three symmetry equivalent con- 
duction bands give rise to a hole pocket at T and electron 
pockets at L (see Fig.[TJa)). Their band dispersions, e v for 
the valence band and (i = 1,2,3) for the conduction 
bands, have been approximated by a parabolic form which 
describes them well near their extrema, in agreement with 
ARPES experiment [9J. 

Our Hamiltonian is composed of a one-electron part Ho, 
containing the one-electron energies, and a Coulomb inter- 
action part W, which represents the electron-hole interac- 
tion between the valence and the conduction bands. Below 
T c , this interaction allows electron- hole pairs, i.e. excitons, 
to form a condensate described by the order parameter A. 
The calculated Green's functions G v and G c describing the 
bandstructure at T and L respectively have the following 
forms [T2] 

G v (p,z)= \ -T\(z-6l(p + Wi)), 
D(p, z) r~ 

-|A| 2 £ ^jmKz-eTip + wj))) 

m,j^i 

(sijm is the permutation symbol). The vectors Wi, called 
spanning vectors, correspond to the distance between T and 
the L points. The denominator D, common to both Green's 
functions, is 

£>(p, z) = (z- e v (p}) Y[(z - e l c (p + wi)) 

% 

-\A\ 2 J2U(z-4(P + ^))- 

i j^i 
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Fig. 1. (a) On the left, schematic picture of the Fermi surface (in 
the plane perpendicular to k z ) considered in the model, placed on 
the Brillouin zone of lT-TiSe2- On the right, side view of the Bril- 
louin zone showing high symmetry points. Graphs (b) and (c) show 
band dispersions with their spectral weight (photoemission spectra) 
calculated with the spectral function at T and L along the high 
symmetry direction FM and AL (long axis of the electron pockets) 
for an order parameter A = meV and A = 100 meV respectively. 
The spectral weight carried by the band is indicated in grayscale. 
In Graphs (c), the thin dashed lines indicate the position of bands 
having negligibly small spectral weight. 

The zeroes of this denominator give the renormalized 
band dispersions, which depends on the order parameter 
A. They are common to T and L, i.e., for G v and G c , 
as one expects for a CDW characterized in our case by 
the spanning vectors wi. The spectral function, A(p, ft) = 
— Im[G(p, ft + i5)]/ir, describes the one-electron spectrum, 
essential for our purposes. It provides us with the spectral 
weight (SW) carried by the dispersions in the process of 
photoemission. Fig. QJb) and (c) present calculated pho- 
toemission spectra at T and L along the high symmetry 
direction TM and AL (long axis of the electron pockets) 
for an order parameter A = meV (normal phase) and 
A = 100 meV (CDW phase) respectively. The SW of 
the dispersions is indicated in grayscale. On Fig. [1] (c), 
the bands indicated by the thin dashed lines have a neg- 



ligibly small SW, so that they do not appear on these 
grayscale graphs. We immediately see that the similarity 
of the bands at T and L is only approximate. At T, as the 
order parameter increases, the lower part of the valence 
band v shifts to higher binding energies, while its top d 
shifts above Ep, opening thereby a gap. Some SW is then 
tranferred from v to c' . At L, with an increasing order pa- 
rameter, the conduction band c does not move but looses 
SW in favor of the backfolded valence band v and the new 
band d (which are the same as those appearing at Y). 




50 100 150 200 250 300 
Temperature (K) 







c,c 






at L or r 



50 100 150 200 250 300 
C) Temperature (K) 





300K 












I-0K- 



-0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 
E-E F (eV) 



Fig. 2. (a) BCS-like temperature dependence chosen for the order 
parameter, (b) Position of the bands at T and L as a function of a 
temperature (for the above order parameter), (c) Evolution of the 
broadened spectra (see text) at L under the effect of temperature. 

In order to introduce temperature effects in the model 
in a simple way, we now choose a particular form for the 

order parameter A(T) = Ao • \J 1 — (T/T c ) 2 where Ao is 
its value at T = OK. This function is drawn in Fig.[2ta) for 
A = 100 meV and T c = 200K. It is similar to a BCS order 
parameter, displaying a steep decrease at T c and a satu- 
ration for T — > OK. Introducing this order parameter into 
the denominator of V and calculating its zeros provides us 
with the curves of Fig. [2] (b). They show the temperature 
dependence of the band positions at T and L. Below the 



transition temperature, the backfolded valence band v and 
conduction band c' shift away from their normal phase po- 
sitions, exhibiting a behaviour very similar to the shape 
of the order parameter. At the same time, the conduction 
band c and a symmetry equivalent version c" stay at their 
inital position. Considering also SWs and broadening the 
J-like peaks with a finite width of 30 meV (for presentation 
purposes), Fig.[2](c) displays over a wide temperature range 
spectra at L (not equivalent to V due to the SWs), where 
the excitonic effects are the most spectacular. Here, below 
Ep, one sees that the evolution of the backfolded valence 
band v is characteristic of the temperature dependence of 
the order parameter (as the backfolded conduction band c', 
which is however not accessible to photoemission, since the 
states are unoccupied). These calculated spectra highlight 
how the real (experimental) temperature behaviour of the 
order parameter can be extracted from ARPES data. One 
sees that the situation is particularly delicate when the or- 
der parameter is small, since the peak of the backfolded va- 
lence band v merges with the peak of the conduction band 
(see arrow on Fig.[2](c)). 

From the condition of conservation of occupied electronic 
states (weighted with the SW), one can also compute the 
temperature dependence of the chemical potential \i. Tech- 
nically, we only take into account states down to —1.0 eV, 
due to the parabolic approximation of the bands around 
their extrema. Then, we compute the number of occupied 
electronic states n occ at T = 300K (taking into account 
the SW of the dispersions) for getting a reference value. 
Then, decreasing the temperature gradually (and possibly 
increasing the order parameter), we adapt the chemical po- 
tential so that the number of occupied electronic states re- 
mains constant n occ (T < 300K) = n occ (T = 300K). This 
has been done for two differents cases. First we fix Ao = 
meV, which means that we look at a normal system, ex- 
hibiting no transition. Then we fix Ao = 100 meV, as be- 
fore, for the excitonic insulator system. Fig. [3] (a) shows the 
behavior of the chemical potential. One sees that without 
excitonic effects, a system having the modelized configu- 
ration of 1 T-TiSe2 undergoes already a chemical potential 
shift of about 12 meV over the range of 300K due to the 
change in the thermal occupation of electronic states. If an 
excitonic phase transition with Ao = 100 meV sets in, a 
drastic change happens around T c and the chemical poten- 
tial shift increases up to nearly 60 meV. Such a large effect 
should be visible in ARPES, but it must be emphasized 
that this holds only for the simplified 1 T-TiSe2electronic 
bandstructure of this model (in particular without consid- 
ering the other Se4p-derived valence bands). 

Knowing the band dispersions around Ep and their shift 
due to the combined effect of the chemical potential and 
the order parameter, it is also possible to estimate the free 
carrier density n = nhoie+ ^electron in the system, composed 
of holes in the valence band and electrons in the conduction 
bands. In Fig.[3](b) (left), in comparison with the normal 
metal (Ao = meV), the excitonic insulator (Ao = 100 
meV) displays a strong free hole density decrease below 
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Fig. 3. Compared behaviours for the cases of a simple semi-metallic 
system (i.e. Ao = meV, in blue) and an excitonic insulator (i.e. 
Ao = 100 meV, in red), (a) Temperature dependence of the chem- 
ical potential, (b) Temperature depence of the free hole (left) and 
electron (right) density, (c) Left graph: inverse free carrier density 
approximating the resistivity (see text for explanation concerning 
the two regimes). Right graph: the measured inplane resistivity of 
!T-TiSe 2 |3] 



T c . This can be explained with Fig. [T] (b) and (c). As the 
order parameter increases from a zero value, a gap opens 
at T and p remains in the conduction band so that the hole 
contribution to conduction vanishes. In Fig. [3] (b) (right), 
the free electron density of the excitonic insulator displays 
also a strong decrease below T c . Indeed, at L, the bottom of 
the conduction band c looses SW, reducing the free electron 
carrier density. In parallel, the chemical potential moves 
upward in the conduction band (this shift is nonetheless 
smaller than the gap produced by the order parameter), 
making available new states in the conduction band c with 
higher SW than those at its bottom. However, this is not 
sufficient to counterbalance the previous effect. 

In transport measurements, the CDW has a strong sig- 
nature, raising the resistivity p below T c as in an insulator. 
Having now an approximate behaviour of the free charge 
carrier density as a function of temperature, we can also es- 
timate the resistivity in the Drude theory by p = m/ne 2 T 1 
with r the relaxation time and m the effective mass. Fig. 
[3] (c) shows 1/n, one ingredient of p, again for Aq = 
meV and Ao = 100 meV. In comparison to the normal 
system (Ao = meV) which displays a nearly constant 
1/n, the excitonic insulator (Aq = 100 meV) exhibits a 



strong increase of 1/n below T c , as expected. Compared 
to the measured inplane resistivity shown in Fig. [3] (c) 
(right), one sees a relatively good qualitative agreement 
above 170K (indicated as the excitonic insulator regime 
in Fig. [3] (c)). At low temperature, our crude estimation 
of the resistivity, p oc 1/n, does not take into account the 
scattering by the lattice represented by the relaxation time 
r in the Drude theory. Indeed, for low temperatures, r 
increases as a power of T, compensating thereby the de- 
creasing n. This generally happens well below the Debye 
temperature, evaluated as Qd = 200K [16J (indicated as 
the lattice scattering regime in Fig.[3](c)). 

Of course, the order parameter of the excitonic insulator 
model could be computed directly in the framework of the 
model, since it obeys to a gap equation similar to that 
of BCS theory. The chemical potential would be obtained 
as well in a self-consistent way. Notwithstanding, such a 
computation is demanding and goes beyond the discussion 
of the present study. 

3. Conclusions 

The temperature dependence of the excitonic insulator 
phase model adapted to lT-TiSe2 [9 12J has been studied. 
This was done in a simple way by assuming a given tem- 
perature dependent order parameter. Then, photoemission 
spectra were calculated as a function of temperature, with 
the help of the spectral function derived from our model. 
They inform us about the signature of the order parame- 
ter in photoemission and indicate how its real (experimen- 
tal) behaviour can be extracted from ARPES experiments. 
Moreover the upturn in the measured resistivity of 1T- 
TiSe2 was explained in terms of the excitonic transition. 
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